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MATHEMATICAL EXISTENCE RESULTS 
FOR THE DOI-EDWARDS POLYMER MODEL 

LAURENT CHUPIN* 


Abstract. In this paper, we present some mathematical results on the Doi-Edwards model 
describing the dynamics of flexible polymers in melts and concentrated solutions. This model, devel¬ 
oped in the late 1970s, has been used and tested extensively in modeling and simulation of polymer 
flows. From a mathematical point of view, the Doi-Edwards model consists in a strong coupling 
between the Navier-Stokes equations and a highly nonlinear constitutive law. 

The aim of this article is to provide a rigorous proof of the well-posedness of the Doi-Edwards model, 
namely it has a unique regular solution. We also prove, which is generally much more difficult for 
flows of viscoelastic type, that the solution is global in time in the two dimensional case, without 
any restriction on the smallness of the data. 

Key words. Polymer, Viscoelastic flow, Doi-Edwards model, Navier-Stokes equations, global 
existence result. 
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1. Introduction. Numerous models exist for describing fluids with complex rhe¬ 
ological properties. They generally are of great scientific interest and have a rich phe¬ 
nomenology. Their mathematical description remains challenging. We are interested 
in this article to a model - the Doi-Edwards model - which was one of the foundation 
of the most recent physical theories but for which the mathematical theory remains 
very poor. 

M. Doi and S.F. Edwards wrote a series of papers [101 HH lEl E] expanding the 
concept of reptation introduced by P.G. de Gennes in 1971. This approach was then 
taken up in a famous book in polymer physics in 1988, see [14]. Since this model 
was derived, numerous studies have been carried out either from a physical point of 
view, either from a numerical point of view. Moreover, several other models were 
born: simplified models using for instance the Independent Alignment approximation 
or the Currie approximation [7], more complex models like the pom-pom model |28] . 
Finally, much progress has been made on the modeling of both linear and branched 
polymers. However, from a mathematical point of view, it seems that no justification 
was given even for the pioneering model. 

Nonetheless we can cite some recent theoretical papers on this subject, see [Hdl], in 
which the authors are only interested in specific cases: one dimensional shear flows 
under the independent alignment assumption in |19] . flows for which the coupling 
between the velocity and the stress is not taking into account, see [4]. More generally, 
there seems to be a real challenge to obtain global existence in time for models of 
polymers. The most caricatural example is the Oldroyd model for which the question 
of global existence in dimension 2 remains an open question, see [24] for a partial an¬ 
swer. However, there exists polymer models for which such results are proved. Thus, 
for the FENE type models, N. Masmoudi [27] proved a global existence result in di¬ 
mension 2. Similarly, for integral fluid of type K-BKZ such results hold too (see m)- 

The aim of this paper is to prove relevant mathematical results on this relevant phys¬ 
ical problem. The first one is the following (a more specific version of this result is 
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given by Theorems 13.11 and 13.21 page [5]) : 

Theorem 1.1. There exists a time t* > 0 such that the Doi-Edwards model 
admits a unique strong solution on the interval time [0,t*]. 

By the expression ’’strong solution” we mean a sufficiently smooth solution so that 
each term of the system is well defined, as well as the initial conditions (corresponding 
to the time t = 0). The lifetime of the solution, i.e. the value of time t*, is not 
easily quantifiable. In practice it is well known that for the Navier-Stokes equations - 
modeling newtonian behavior, the question of long time existence is still an open one. 
For the Newtonian fluids, the only existence results, for long time and for any data, 
correspond to the two dimensional case. However, we have above pointed out the 
difficulty to get this kind of result even in 2E for some viscoelastic fluids. The major 
point of this paper is the proof indicating that the model of Doi-Edwards admits a 
strong solution for long time in 2D (a more specific version of this result is given by 
Theorem 13.31 page [9]): 

Theorem 1.2. For all time t* > 0, the Doi-Edwards model admits a unique 
strong solution on the interval time [0,t*]. 

The paper is organized as follows. First - in Section [U we introduce the Doi-Edwards 
model specifying the physical meaning of each contribution. This second section 
ends by a dimensionless procedure that allows us to write the model with only three 
parameters (the Reynolds number, the Weissenberg number and the ratio between 
solvent viscosity and elastic viscosity). In Section |31 we present the mathematical 
framework as well as the assumptions which are physically discussed. The main results 
are given at the end of this section. Section |4] is devoted to fundamental preliminaries 
which correspond to some key points of the next proofs. The first two provide a 
priori bounds which will imply that that the stress defined in the Doi-Edwards model 
is automatically bounded. The third preliminary give a Gronwall lemma with two 
time variables. The fourth preliminary result is about the maximum principle which 
can be applied many times to estimate the memory of the fluid. The last preliminary 
result is about a Cauchy problem arising in the global existence proof. The three last 
sections dSHH and [71) are devoted to the proof of the three mains results , namely: the 
local existence result in Section the uniqueness result in Section [5] and the global 
existence result in Section [71 Some open questions are presented by way of conclusion. 

2. Governing equations. 

2.1. Conservation laws. In this paper we are interested in the flow of isother¬ 
mal and incompressible fluids. The incompressibility implies that the mass conser¬ 
vation is equivalent to the free-divergence of the velocity field. The isothermal as¬ 
sumption implies that only one other conservation law suffices to describe the flow: 
it corresponds to the law of conservation of the momentum (Newton’s second law of 
motion). This equation is written as a balance between the material derivative of the 
velocity and the divergence of the Cauchy stress tensor. For a polymeric liquid, the 
equations of conservation can hence written as a system coupling the velocity field v, 
the pressure p and the extra-stress tensor cr: 


p dtV -1 Vp — = div cr, 

diw = 0, 
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where p is the fluid density and rjs the solvent viscosity. The notation dt corresponds 
to the material derivative dt = dt + v ■ V. 

2.2. Constitutive equation. A fundamental result of the Doi-Edwards theory 
is the expression for the stress tensor cr which applies when the chains which compose 
the fluid are relaxed within their tubes. More precisely, the stress can be deduce from 
a tensor S denoted the orientation order parameter of the chains. Although the chain 
tension is permanently at the equilibrium value, the orientations become anisotropi- 
cally distributed as a consequence of the flow, and a stress develops accordingly. The 
stress is then modelized by (see P page 2056]): 

cr{t,x) = -^ J ^ S{t,x,s)ds, 

where Ge is a characteristic modulus and i is the equilibrium value of the contour 
length of the chains. The quantity s is a arc-length coordinate along the primitive 
chain. 

To evaluate this tensor S, M. Doi and S.F. Edwards write S = ^d). Here, u is 

a unit vector along the tangent to the primitive chain which depends on time t, spatial 
position X and length s, and S denotes the identity tensor. The entire d corresponds 
to the dimension of the spatial coordinates (in practice d = 2 or d = 3). The average 
is over the distribution of these vectors in the ensemble of chains, i.e., more explicitly 

S{t,x,s) = r f{u\t,x,s)[u<S>u—^S)du, 

Jgd-i d 

with /(u; t, X, s) given the orientation distribution function. To obtain an expression 
for this distribution, the history of motion must be found. To this purpose, let us first 
recall the relevant aspects of the Doi-Edwards model. They are 

■/ The polymer moves randomly inside the tube executing one-dimensional 
Brownian motion. Moreover tube segments are randomly oriented when they 
are created and deform affinely thereafter; 

/ If the system is macroscopically deformed, the polymer conformation is also 
changed as presented on Figure [2Tl 

■/ The macroscopic motion and Brownian motion coexist, independently from 
one another. 



V 


Fig. 2.1. (inspired from the book JP) page 2058]) - When a macroscopic deformation is applied, 
a polymer chain is transformed into a new conformation. The new chain is on the curve which 
is the affine transformation of the initial curve. The new position of each segment is obtained by 
retraction, preserving the initial lengths. 
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All these considerations being taken into account, it is possible to write (see m page 
277]) 

^ + 00 

S{t, x,s) = - dTK{t, T, X, s) T, x)) dT. (2.1) 

Jo 

It makes appear the deformation gradient tensor G which depends not only on the 
current time t and spatial variable x but also on an other time T. The time T allows 
to take into account all the history of the motion. The deformation gradient G fulfills 
the differential equation (see my- 

dtG + drG = GVv. 

Finally, the integral kernel K satisfies (the coefficient is a curvilinear diffusion 
coefficient): 

dtK + drK + (Vv: ^ S ds^dsK - Ded'^.K = 0. (2.2) 


Remark 2.1. 

1. Note that the equation (EB) does not explicitly defined the orientation tensor 
sinee S again appears in the equation (113). It is one of the difficulties to 
obtain existence results. 

2. The usual formulation use the time t and a other time t' in the past. Morally, 
the deformation gradient tensor G measures the deformation between these 
two times. In the present paper we select as independent variable the age 
T = t — t', which is measured relative to the current time t. This viewpoint 
is relatively classical in the numerical framework for integral models, see for 
instance [WA [HJ 

The model is closed with the expression of the function see [H eq. (7.115)]: 


S^{G) 


1 /{G ■ u)®{G ■ u)\ 1 

<1G-m1>o\ \G-u\ /o d 


(2.3) 


where the brackets (•)o correspond to the average over the isotropic distribution of 
unit vectors u e 


2.3. Boundary and initial conditions. The previous equations are supple¬ 
mented by boundary and initial conditions. Throughout this article we restrict to the 
case where the macroscopic field is assumed to be periodic. Thus the only condition 
that we impose on the unknowns v, p, K and G with respect to the variable x is to 
be periodic. Clearly this “simplification” is purely mathematical and it will be inter¬ 
esting to treat a more physical case imposing, for instance, the value of the velocity 
at the macroscopic boundary. 

By definition of the integral kernel K, we impose the following conditions: 

£ £ 

K{t,0,x,s) = l and K{t,T,x,--) = K{t,T,x,-) = 0. (2.4) 

In the same way, the quantity G(t, T, x) which corresponds to the deformation gra¬ 
dient from a past times t — T to the current time t must naturally satisfies 


G{t, 0, x) = 6. 


(2.5) 
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The initial conditions correspond to the value that we impose at time t = 0. We 
assume that we know the velocity at this initial time and at any point x of the 
domain. In the same way we assume that we know all the past of the flow before the 
initial time: we then know the value of G and K at t = 0, for any age T and at any 
point X. To summarize, we assume that there exists an initial velocity Vq, an initial 
deformation Go and an initial function Kq such that 

v{0,x) = Vo{x), 

G{0,T,x) = Go{T,x), (2.6) 

K{0, T, X, s) = Ko{T, X, s). 


2.4. Remark: The I.A. approximation. A common approximation for the 
Doi-Edwards model, called the independent alignment approximation (LA. approxi¬ 
mation, see m section 7.7.2]), is to neglect the transport term (Vv : S ds^dgK in 
the equation (12.21) . and also to simplify the expression (12.31) of the function ^ using 

cAiA) _ / (G ■ m) (X) (G ■ m) \ _1 

^ ^ \ |G-m|2 a d 

The LA. approximation is actually quite popular in the rheology literature (see e.g., 
[2l[23l[26]) since the corresponding configurational equation for K can be explicitly 
solved using the Fourier series. The stress tensor cr is then more simply given by 
(see [m Equation 7.195]): 


(T^^^\t,x) = Ge 


r-l-GO 

Jo 


m{T)y^^^'>{Git,T,x))dT, 


where 


m{T) = 


s 

p odd 


8£>e 

TT^ P 


exp 


-TDe 

p 



(2.7) 


For such a model the global existence result is a consequence of a general result on 
viscoelastic flows with memory, see [3]. Nevertheless, it is also well known that this 
approximation causes serious error in certain situations, this is clearly specified in the 
seminal book M- More precisely, it is proved that LA. predicts a negative Weis- 
senberg effect (see [18]) while the version without LA. predicts a positive Weissenberg 
effect (see [25]). To paraphrase M. Doi [9] page 2064]: ’’Mathematically [...] there 
seems no a priori reason why the term (Vv : ds)ds can be neglected compared 

with the term Dgd^”. 

2.5. Dimensionless procedure. In order to recover characteristic properties 
of the system, we use a nondimensionalization procedure. We denote by L a char¬ 
acteristic macroscopic length, by E a characteristic velocity of the flow. It is then 
natural to define a dimensionless coordinates x*, a dimensionless velocity v* and a 
dimensionless time t* by the following relations 

X = La;*, V = Vv*, t = y'^* ■ 

For polymer flow, there exists also two microscopic characteristic sizes which corre¬ 
spond to the length I and to the diffusion coefficient Dq. They allow to define a 
dimensionless microscopic length s* and another dimensionless time T*: 

P 

s = ls\ T=—T\ 
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Finally, in a dilute polymer solution, two viscosities naturally appear: the solvent 
viscosity rjs and the elastic one defined using the characteristic modulus: r]e = 

If we denote hy r] = % + r]e the total viscosity, then we defined the dimensionless 
pressure and stress as follows 



cr = - (T . 

L 


Taking into account all these new unknowns and new variables, the complete system 
reads (without * in the notations): 


fHe dtv + Vp — (1 — uj)Av = div cr, 


(2.8a) 


diw = 0, 

cr(t, x) = u! 


1 

J ^ S{t,x,s) 


ds. 


r+OD 

S{t, X,s) = - dTK{t, T, a;, s) T, x)) dT, 

Jo 

dtG+ ^dTG = G-Vv, 

2I(e 

<^tK+^^dTK+ (v-u: J 




(2.8b) 

(2.8c) 

(2.8d) 

(2.8e) 

(2.8f) 


In this set of equations, Die is the usual Reynolds number, w stands for the viscosities 
ratio and SUe is the Weissenberg number defined by the ratio between the macroscopic 
time and the microscopic time. More precisely we have 


SHe = 


pVL 


U) = 


Ve 


me = 


eiD, 


Tj rj LjV 

The functions 5^ is always defined by the relation (12.31) . 


The goal of the rest of the paper is to analyze, from a mathematical point of view, 
the existence of a solution to the system (12.81) . More exactly, by given initial data 
(dq, Go; is there a triplet of functions (v, G, K) which coincides with the data at 
initial time and such that the previous system holds for any future time? 

3. Mathematical framework, assumptions and main results. 

3.1. Notations. The integer d stands for the spatial dimension of the flow. It 
will be equal to 2 or 3 in the first parts and exclusively equal to 2 in the Section [7] 
where we prove a global existence result. 

Notations for functional analysis - 

— The d dimensional torus is denoted T. 

— For all real n > 0 and all integer g > 1, the set W"’'^ corresponds to the usual 
Sobolev spaces with respect to the space variable a; e T. We classically denote 
L® = and R" = and we do not take into account the dimension 

in the notations, for instance the space (IF^’'^)^ will be denoted 
— All the norms will be denoted by index, like ||u||^yi,<,. 

— Since we are interested in the incompressible flows, we introduce 

Hg = {v e L^.; divu = 0}. 
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— The Stokes operator is introduced, with domain S>{Aq) = n Hg, 

whereas we denote (see [SI Section 2.3] for some details on this space) 

f+oo 

= {ve Hq; \\v\\^g + ( \\Aqe~*^'>v\\l, dt) < +oo}. 

Jo 

— The notations of kind L’'(0,t*; X) denote the space of r-integrable functions 
on (0, t*) with values in the space X. For instance G e T'’(0, t*; L”L|) means 
that 

l|G'llL''(ot*-L'=°L'’) ■= f sup ^ r \G{t,x,T)\’^dx\ dt<+co. 

Notations for tensorial analysis - In System (12.81) . the first equation (I2.8al) is 
a vectorial equation (the velocity d is a function with values in R*^). The equa¬ 
tions (I2.8cl) . (I2.8dl) and (I2.8el) are tensorial equations (the stress cr, the orientation S 
and the deformation tensor G are functions with values in the set of the 2-tensors). In 
the following proofs, we need to work with the gradient of such 2 -tensors, that is with 
3-tensors, and even with 4-tensors. We introduce here some notations for tensors. 

— The set of linear applications on the d-dimensional space is denoted £(M'^). 
— The products A®B^ A - B and A : B between two tensors of order p and q 
are respectively defined component by component by 

(A^B). . = Qi, i bi, i , 

where we use the Einstein convention for the summations with indexes k 
and i. 

— Note also that all these products are inner products on the set of the p-tensors. 
It allows us to define a generalized Froebenius norm: 

1^1" = S °'l,...,^p■ 

if,... ,ip 

We conclude this section introducing the constant C. This constant stands for any 
constant depending on the data of the problem: initial conditions, physical parame¬ 
ters... In some cases, informations will be given on the dependence of this constant 
(see for example Section |3 where we explain that this constant may depend on time t* 
but must remain bounded when t* is bounded). 

3.2. Assumptions. The results proved in this article requires some assumptions 
about the data. In addition to the assumptions on the regularity of the initial condi¬ 
tions that will be specified in each theorem statement, we will need some ’’natural” 
assumptions. 

/ The first assumption relates to the initial deformation Gq: 

Sy > 0 ; detGo ^ 7 . (3-1) 

We note that in many applications, the fluid is assumed to be initially quiescent. In 
that case, we have Gq = S and detGo = 1- Moreover, we will see (equation (14.Sp 
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in the preliminary section that the quantity detG is only convected by the flow. 
If the fluid is assumed to be at rest in the past (that is for T large enough), then 
we always have detGo = 1. The assumption on the positiveness of det Gq allows us 
consider, for instance, such cases. 

/ The second assumption relates to the initial memory mo = —BtKq: 


mo ^ 0, (3.2a) 

drmo ^ 0. (3.2b) 

The assumption ()3.2a|l corresponds to the fact that the quantity mo describes the 
memory of the fluid, that is the weight that must have the quantity =5^ in the flow 
via the relation (I2.8dp . It is physically positive. The assumption (13.2b)) indicates that 
the memory decreases with the age T: It is linked to the principle of fading memory, 
see [5]. 

In the integral models, that is to say when the memory is explicitly given in terms of 
age T, it is a combination of exponentially decreasing functions (see for instance the 
Doi-Edwards model under the LA. approximation, subsection 12.41 and more precisely 
the expression (EZD of the memory). Such decreasing behaviors will be prescribed 
in the functional spaces with exponential weight. For example, we will impose that 
there exists /i > 0 such that mo ^ C e“™' . 

3.3. Main results. The first result concerns an existence result for strong solu¬ 
tion. It is a local in time result: 

Theorem 3.1 (local existence). Let r e]2,+co[, q e]d,+co[ and fj, > 0. 

If the data Vq, Gq and Kq satisfy the assumptions (13.11) . (I3.2al) and have the following 
regularity 

Vo e Go e BtGo e L^Ll, 

then there exists t* > 0 and a strong solution (u,G,K) to System (12.81) in [0,t*], 
which satisfies the initial and boundary conditions (1^ . (1^ and (1^ . Moreover 
we have 

VG L’'(0,t*;IEj’'?), 

GGL“(0,t*,L®IEi’‘?), 

K, drK, e'^ drK g L”(0, t*; L®L® J, 

dTViLGL®(0,t*;L^L?_,nL|L«^J, 

and the memory m = —drAT remains non negative. 

Remark 3.1. In this article, we will not give any result on the pressure p. In 
practice, the latter is regarded as a Lagrange multiplier associated to the divergence 
free constraint. It can he solved using the Riesz transforms. More precisely, taking the 
divergence of the first equation p.8a|l of System (12.81) we use the periodic boundary 
conditions to have 


dtveM{QX-Ll), 
dsG, dtGe L^{0,t*;Lif;Ll), 
dtK e L^{0,t*-L^Ll,), 
d,KeL^iO,t*-,L^Ll,), 


p=—{—A) ^ div div{cT — V ®v). (3.3) 

From Theorem no we can prove that the solutions of System (12.81) discussed in this 
paper have cr~v®v in L®(0,<*;L^). The pressure in the solution of (12.81) is meant 
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to be given by (lO) . 

We will see during the proof of this theorem 13.11 that one of the key point is the 
behavior of the memory m = —dxK for large value of the age T: if the memory 
is exponentially decreasing at t = 0 (with respect to the age variable T) then the 
solution will be exponentially decreasing for any time t > 0. 

In the same way, it is possible to prove that if the memory m is initially decreas- 
inf3 (with respect to the age variable T) then the solution will be decreasing for any 
time t > 0. The proof - which is not given in the proof of Theorem 13.11 - consists in 
derivating twice the equation (l2.8fD with respect to T, and next in applying the max¬ 
imum principle (see the subsection 14.41 page ITS)) to the function cItto. Therefore, the 
assumption (I3.2bl) . which is not necessary to obtain local existence, is also preserved 
in time. 

We will show that the solution obtained in Theorem 13 .1 1 with this additional assump¬ 
tion (|3.2bl) is the only one in the class of regular solutions. Precisely, the result reads 
as follows. 

Theorem 3.2 (uniqueness). Let t* > 0. 

Let (ui,Gi, Ki) and {u 2 ,G 2 , K 2 ) be two solutions to System (12.81) satisfying the 
initial and boundary conditions (EH), (E^) and (1^ . If we have, for i e {1,2}, 

V«,eL2(0,l*;L”), 

G,eL^{0,t*-L^{L^nWy)), 

drK, e L“(0,r;L?L”, n 

and if each rui = —drUi is decreasing with respect to T then the two solutions coincide. 

Obviously, the solution obtained in Theorem 13.11 satisfies the regularity requested 
in (13.4p . Combining Theorems 13.11 and 13.21 we get a local and uniqueness result. 
In the two-dimensional case, it is possible to show that the solution {v,G,K) of 
problem (12.811 exists for any time t* > 0. More precisely we have the following result: 

Theorem 3.3 (global existence in 2D). Let r e]2, -l-oo[, q e]2, -l-oo[ and p, > 0. 
We assume that y + “ < ^ that the data Vq, Gq and Kq satisfy the same as¬ 
sumptions that in Theorem \8.1\ 

Let t* > 0 be arbitrary. 

There exists a constant C depending only on the data with C bounded for bounded t*, 
and a solution {v,G,K) of (12.81) satisfying the initial and boundary conditions (12.41) . 
(I2.5[l and (1^ such that 

^ C', 

The estimates (|3.5|1 announced in the theorem |3.3l above are sufficient to prove a global 
in time existence of a solution {v,G,K). In fact if (I3.5|l holds then it is possible to 
prove - principally using the lemmas introduced in the proof of the local existence 
result - that the solution {v, G, K) of (12.81) at time t* have the same regularity that 

^ Be careful not to confuse the terms ” exponentially decreasing” and ” decreasing”. The first 
means that m is bounded by a function of the form e“^, while the second means that dTcn ^ 0. 
Moreover the first is a global property, while the second is a local property. 




(3.5) 
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at time 0. Applying the local result (theorem 13.IL we deduce that the solution can 
not blow up in finite time. 

4. Preliminaries. In this section we give some results which will be using during 
the different proofs of the previous theorems. 

4.1. Some bounds for the function . One of the key points of the proof of 
global existence lies in the fact that the stress cr, defined by (I2.8cl) - (l2.8dl) is bounded. 
The first result in this direction is the following result concerning the function 

Proposition 4.1. The function defined by the relation (12.31) is of class 
on >C(R‘^)\{d} and satisfies the following properties: 

> 0 ; VG e £(R‘^)\{d} , \J^{G)\ ^ (4.1a) 

3.^^ ^ 0 ; VG e , \G\\Y'{G)\ ^ (4.1b) 


Proof Recall the definition (12.31) of the function 

^ 1 / (G-u)®(G-u) 

^ ^ <|G-u|>o\ |G-u| 

•/ We first notice that the following inequality is obvious 




l</>o| = f /(M)dM s: r |/(M)|dM = <|/|>o, 


(4.2) 


so that the first point of the proposition 14.11 is a direct consequence of the inequality 
|A(x)B| ^ |A| |B| for all tensors A and B, and of the relation |(5| = Vd. More precisely 
we obtain, for all G e £(IR‘^)\{Q}: 


\^{ G )\^1 + 


Vd’ 


which corresponds to (I4.1al) . 

/ For the second point (I4.1bl) . we write the tensor ^{G) component by component: 
for any (fc, £) e {1,..., d}^, 


^(G)fc, 


1 /(G ■ u)k (G ■ u)£\ 1 

<|G-m|>o\ |G-m| /o d 


It is then not difficult to evaluate each components of the derivativ^l tensor S^'{G). 
For instance for {i,j,k,£) e d}^ we have 


Sg,A\G-u\) 


(G • u)i Uj 
\G-u\ 


^The application being defined on an open set of with values in its differential 

is an application with values in £(£(IR^),Consequently, ^'{G) can be identify to a 4-order 
tensor whose components are = ^Gij^{G)M' 
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More generally we obtain 

dG,,nG)ke = 


-1 


' {G ■ u)iUj\ / {G ■ u)kiG ■ u)i' 



|G-m| /o\ |G-m| / 

1 1 2 


(\G-u\)o' 

\ \G-u\ 

+ .. < 

/ (G • u)k{G ■ u)i{G ■ u)iUj \ 


Taking the norm in equality (14.3p and using (j4.2L we deduce 

|^'(G)| ^ 2(1 +Vd)- 


l«l>o 


\G-u\}o- 


(4.3) 


(4.4) 


We note that (|m|)o = (l)o = r(d/ 2 ) ■ application 

G 1 -^ (|G • m|)o is a norm on the finite dimensional space £(R'*). This norm is then 
equivalent to G |G|: there exists a constant Cd such that for all G e £(]R‘^) we 
have 


(\G-u\)o^Cd\G\. 


The inequality (14.4p becomes 

^ 27r"*/2 

|G||.5''(G)| « 2(1 + 

that concludes the proof of the Proposition 14.11 □ 

4.2. Positive norm for the deformation G. In practice, we will prove that 
for any solutions {v,G,K) to the System (|2.8I) the deformation gradient tensor G has 
a positive norm. More precisely we have 

Lemma 4.2. Let v he a free divergence vector field on T and G be a solution to 
Equation (I2.8ep with initial conditions satisfy detG|t=o ^ 7 > 0 and G\t=o = d . 
There exists a constant 7 such that for all (t^T,x) e (0,t*) x R+ x T we have 

|G(t,T,(r)| ^ 7 > 0 . 


Proof A simple calculation shows that the quantity det(G) satisfies 

X>det(G) = divv det(G) = 0, (4.5) 

where V refers to the one order derivating operator V = dt + The value det(G) 

is then constant along the characteristic lines. Since all the characteristic lines start 
from the lines {t = 0} or {T = 0} we deduce from the assumptions that det(G) > 
min( 7 ,1) on (0,t*) x R+ x T. 

Due to the inequality of arithmetic and geometric means, we have 

|G|2 = Tr(^G- G) > 2^det(^G- G) = 2|det(G)| > 2 min( 7 ,l). 


That concludes the proof of Lemma 14.21 taking 7 = y/2 min( 7 ,1). 
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For the local result we will use the fact that the derivative S^'{G) is bounded. That is 
clearly not the case if we only use the Proposition l4.ll We then introduce the function 

i^(G) = ^(G)x(|G|), 

where x e [0,1]) satisfies xl[o, 2 ] = 0> xl[ 7 ,+oo[ = 1 and x' < f (see the 

Figure SH]). 



Fig. 4.1. An example of truncation function x- 


Using the Proposition 14.II it suffices to derivate 

y'{G) = y'{G) x(|G|) + no) (X) yli x'(|G|), 

to deduce 


VGe£(R‘') \y’{G)\^yao and '(G)| ^ .5^^, 


where + So5^oo- Using the Lemma |4?^ we have the following consequence: 

Lemma 4.3. Under the assumption (EH), the solution of the System (12.81) is 
the same if we use instead of . In other words, we can assume that S^'{G) is 
bounded on (0,t*) x 1R+ x T. 

4.3. A Gronwall lemma with two times. The choice that was made in this 
article is to use the current time t and another time T corresponding to the age of 
the flow. This choice simplifies the expression of the orientation tensor S since the 
t-dependence does not appear in the integral bounds (see (l2.8dl) L The price to pay 
is that the derivatives in time in the evolution equations of K and G involve both t 
and T. In this framework, the following lemma is the analog of the classical Gronwall 
lemma for functions depending only t. Its proof is based on a change of variable that 
easily allows to follow the characteristics. It is proved in [2]. 

Lemma 4.4. Let f : K"*" i—> M"*" a positive and locally integrahle function. If a 
function y : K"'" x 1R+ i—> M satisfies, for all (t,T) e (0,t*) x ]R+; 

dty{t,T) + < f{t)y{t,T) 

then we have, for all (t,T) e (0,<*) x IR+; 

yit,T) C{t,T)exp(^^ f{t')dt' 


where f{t, T) 


2 /( 7 ^-^, 0 ) ift^WzT, 
y(0,t-2IIeT) ift>WzT. 
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4.4. A maximum principle. The memory of the fluid is described by the 
function K. This function satisfies the equation (l2.8fD . The form of the equation (I2.8f|) 
is particuiar in that it checks a ” maximum principie”. This wiii be used repeatediy 
in the foiiowing sections. 

Lemma 4.5. Ifdgg e T”s) and v is free divergence on T then the solution 

f(t,T,x^s) to the following system 


dtf H- drf + g dsf - dif = 0, 


/L=o = /o’ 


/|t=0 


/I 



(4.6) 


satisfies the following maximum principle on (0,t*) x ]R+ x T x ^): 


min{ inf /o, inf /i} < / < max{ sup /o, sup /i}. 

T.x^s t,x,s T,x,s t,x,s 


Proof Considering / — min{inf /g, inf /i} or max{sup /g, sup /i} — / instead of /, 
it suffices to show that the soiution / is non negative if the data /g and /i are non 
negative. 

We muitipiy the first equation of (14.61) by the negative part f~ of /, and we integrate 
with respect to x and s. It is important to notice that the function / does not 
necessary satisfy the boundary conditions f \ _i = f \ _ i = 0, but that its negative 

part f~ satisfies these conditions. Using integrations by parts, we deduce 


StWf 


-||2 
LI 




+ 5si^.rii= 



Ssg\f P- 


Since dgg e we obtain 


IIlJ . < Wily 

The Gronwaii iemma with two variabies (see the iemma 14.41) impiies that 

Wf~WliS''^^T) ^ C(C'r)exp(||(5s5||Li(0.t;L»j), 


where the function f oniy depends on the vaiues of ||/ ||^2 on the boundaries {t = 0} 
and {T = 0}. In the case where the data /g and fi are non negative we have C = 0 
and we deduce that ||/“|||2 = 0. We conclude that / is non negative too. □ 


4.5. A Cauchy problem involved in the proof of global existence. In 

this section, we are interested in the following Cauchy problem 

|y(0) = 0, y'(0)=ei, 

where ^g, ^i, ^2 and k are positive constants. 

This problem will occur during the proof of global existence theorem, section 0 More 
precisely, we will see that with a good choice of parameters ^g, ^ 1 , ^2 and k, we 
can control for arbitrarily long time. This will allow to control the 
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stress <J. 

Although this equation (14.71) is nonlinear and of order 2, we can explicitly give the 
solution. We will see that its expression makes appear the function F : M+ —> R+ 
defined by 


F(A) = r dx. (4.8) 

Jo 

Clearly this function is one-to-one (increasing) from to [0,£[, where the real i 
denotes the limit: 

i= \ jj;. 

Jo 

Proposition 4.6. The Cauchy problem (14.71) admits a unique solution given by 


y{x) 




e c 1 

'= + 1 x) 


for all 


X G 



- -k + 1 

r+1 


[■ 


(4.9) 


Proof. The equation (ITTl) is a two order ordinary differential equation and it is 
possible to apply the Cauchy-Lipschitz theorem: there exists a unique local solution. 
Moreover, we clearly have y” > 0 and then y' ^ f,i > 0. The equation (ITTI) also write 

Using boundary conditions given in (14.71) we obtain the following first order ordinary 
differential equation: 


/ ^ _ _5o_f„+f„A+i 

y = fc+i 

Making appear the function F, we deduce that 

M ' 

= fie '=+1 . 

Since F(0) = 0 and y(0) = 0 we integrate and deduce 



The expression (gH) given in Proposition 14.61 follows since F is one-to-one. □ 

5. Proof of the local existence result: theorem 13.11 

5.1. Strategy: a point fixed formulation. In order to prove the local ex¬ 
istence result, we rewrite the set of equations (12. 8p as a fixed point system. More 
precisely, we consider the mapping (the spaces will be further given) 

$ : {v,G,K), 


dehned as follows: 








Existence for the Doi-Edwards polymer model 


15 


/ Velocity problem - The velocity v{t, x) is the solution of the following Stokes 
problem 


Dte dtv + Vp — (1 — w) At) = /, 

< divt) = 0, (5.1) 

^ v\t=o = Vo, 

where the source term f contains the nonlinear term of the Navier-Stokes equations 
and the term coupling velocity and stress, namely / = —Diet) • Vv + diver. The 
stress cr{t, x) is defined by 



S{t, X, s) ds, 


where the orientation tensor S{t, x, s) is given by 


(5.2) 


S{t, X, s 



drKit, T, X, s) y{G{t, T, x)) dT. 


(5.3) 


/ Deformation problem - The deformation gradient tensor G{t, T, x) is the 
solution of the linear equation 


dtG + v-VG+ —drG = G • VU, 


G\t=o = Go, 


G|t=o = 


(5.4) 


•/ Memory problem - The scalar quantity K{t,T,x, s) is the solution of the 
following linear equation 

dtK + v-yK+ + gdsK- = 0, 

^ = 0 , ( 5 - 5 ) 

^ K\t=o = Ko, K\t=o = 1 , 

where g, which only depends on t, x and s, is given by 

g{t,x,s) = Vv{t,x) ■. r S{t,x,s')ds'. 

Jo 

The goal of the next subsections is to analyze these problems (15.111 , (15.4|) and (15.5|) 
independently. We will see at the end of this analyze that it is possible to apply the 
Schauder fixed point theorem for the function $ with adapted functional spaces in 
order to deduce the result of Theorem 13.11 

5.2. Estimates for the velocity v solution of a Stokes problem (15.11) . 

The results for the Stokes system m are very numerous. In this subsection we 
only recall, without proof (we can found a proof in m), a well known result for the 
time dependent Stokes problem. In order to simplify expressions, we use the following 
norm on the velocity field: 
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The result on the Stokes problem (EH) states as follows: 

Lemma 5.1. Let t* > Q, r e]l, +oo[ and q e]l, +oo[. 

If vq g and f G then there is a unique solution v G L'^{0,t*; S>{Aq)) 

such that dtv G L’’(0, t*; to System (15.11) . This solution satisfies 

where the function Fi depends on r, q, uj, Die and the initial value Vq. Moreover this 
function is continuous and nondecreasing on K”*". 

In practice, the function Fx may be chosen as (see m) 

Fx{X) = ^^^{^t\\v4^.,,+X). 

i — UJ ® 

5.3. Estimates for the deformation gradient G, solution of (15.41) . The 

existence and regularity for the deformation gradient G is less classical. As previously, 
we introduce an adapted norm, namely for the deformation gradient 

IIIG'lb := + II tic'll Lr(0,i*;L:»L’) + l|^TG||ir(0^t*;L^L9). (5.6) 


Lemma 5.2. Let 0 < t* < 1, r g]1, +cx)[ and q e]d, +oo[. 

If Go G LflWf’fi drGo e and v g L®(0, t*; ITi'?) n 1.1(0, t*; ITj’?) is free 

divergence then the problem (15.41) admits a unique solution G G L'^{0,t* ; 
such that dfG and BtG belongs to L''(0, t*; This solution satisfies 



where the function F 2 depends on r, q, SUe and the initial value Go- Moreover this 
function is continuous and nondecreasing on R+. 

The proof of a very similar result is given in [2]. One of the differences is that we 
show here dtG G L^{0,t*; L^LD and not only dtG G L’'(0,t*;L®) for any T G R+. 
This difference takes its significance when we will give sense to initial conditions, see 
the remark [01 

Proof The existence of a unique solution to (15.4p follows from the application of 
the method of characteristics (see [El Appendix p. 26]). In practice, the following 
estimates will be made on regular solution G„ which approaches the solution G 
when a regular velocity field Vn approaches the velocity v. The regularity of these 
solutions Gn with respect to t and T comes from the Cauchy-Lipschitz theorem. For 
sake of simplicity, we omit the indexes ” n”. In the following proof, we refer to [15] for 
the passage to the limit n ^ + 00 . The rest of the proof of Lemma lA^ is split into three 
parts: in the first one (see the subsection l5.3.II) we obtain a first estimate concerning 
the regularity of G, and in the subsection 15. 3. 31 we obtain the estimate for dtG. This 
estimate requires an estimate for BtG, which is given in the subsection 15.3.21 

5.3.1. Estimate for the deformation gradient G. Let q > d. We take the 
inner product of the equation (15.41) by 9 |G|'^“^G, and next we integrate for a; G T. 
Due to the incompressible condition divl7 = 0, we obtain 

^g||VU||L»||G||«,. 
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Then we use the continuous injection ^ holds for g > d and making appear 
a constant Cg- 

St\\G\\l, + < gC«||Vu||^i„||Gr^,. (5.7) 

Now, we take the spatial gradient in (15.41) and compute the inner product of both 
sides of the resulting equation with g|VG|‘^“^VG (we will note that this is a inner 
product on the 3-tensor, defined hy A :: B = ai^^k After integrating for a: e T 

we obtain 

dt||VG||«, + ^dT||VG||«, ^ iVGnVul + \G\\VG\‘>-^\y^v\. 

Using the Holder inequality and the continuous injection again, we deduce 

^ 3ga||Vu||^n,||G||^^„. (5.8) 

Adding this estimate (15.81) with the estimate (15.7L we obtain 

^t\\G\\yyl,q + ^^^T||G||,^n<J ^ AC s\'^ v\y^l,q\G\^l,q . 

Using the initial conditions we have 

ll^llw,J’®L=o ~ ll^llwi’'* It=o ~ 

so that the Gronwall type lemma |T4] (see the preliminary subsection l4.3|) implies that 
for (t, T) e (0, t*) x R+ we have 

||G||^c<,(t,T) ^ C(i,7^)exp(4G, II , (5.9) 

Vd 

The assumption Gq e implies C, e L“(0,t*; L®) with 

IICllL‘»(0,t*;L^) ^ mSX I ||Go , Vd}. 

The relation (15.91) now reads 

ll^llL<»( 0 .t*;L:»Wi'‘') ^ IICIlL'=°( 0 ,t*;L”) exp(4Gs||(5.10) 

5.3.2. Estimate for the age derivate drG. We first remark that the deriva¬ 
tive G' = BtG exactly satisfies the same PDE that G (see the equation of (|5.4I) : that 
is due to the fact that v does not depend on the variable T). We then deduce the 
same kind of estimate that dSH): 

MG'Wli + ^^tIIG'IUs < G,||VUi|^yi,<,||G'|U,. 

But the initial conditions differ as follows: 

G'|t=o = 5TGo and G'|t=o = QUe Vu. 


where C(t, T) = 
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This last condition is obtained using T = 0 in the equation (15.41) . Note that this 
result is valid because we are working on regular solutions G„ (see the introduction 
of this proof) such that dtGn is continuous at T = 0. From Lemma [4.41 given in the 
subsection 14.31 we obtain for all {t,T) e (0,t*) xK+ the estimate 


\\G'\\Llit,T) ^ C'(i,T)exp(a JjlVUll^c,), 


(5.11) 


where C{t,T) = | 
For each t e (0, t*) 


||<5rGo|L.(T-^) ift^WcT, 

WelVvllLg(t-WeT) if t>mT. 

we estimate the L®-norm of the function T 


as follows: 


||C'l|L”(t) sS max(2IIe || Vu||L<»(o,t;L’), \\StGo\\l^li)- 


Taking the norm in L’'(0, t*) we deduce that e L^(0, t*\ L”) with 

IIC^I|L'-(o,t*;L“) < max(2Ile || Vu||i<x>(o,t*;L9), ||(?tGo||l:»l« ) t*. 

By assumption we have t* < 1 so that < 1. The relation (jS.llI) now reads 


l|G'||Lr(o,t*;L®L«) < I|C1 |l'-( 0,exp(Gs||Vu||^i(p_j..^CPj). (5.12) 

5.3.3. Estimate for the time derivate dtG. Isolating the term dtG in the 
equation (15.41) we have 


\\StG\\Li 


^ T^WG'hl + N|l-||VG|U. + ||G|U^||VU||i. 

UUe 

^ 2Gs||u||^^c9||G||^c<j. 


Taking the L’'(0, t*; L”)-norm for the variable {t,T), we obtain 


Ild.GI 






^\\G'\\Lr{0,t*-,L^Li) 

+ 2Gs Ill’ll II G||^ai(o,t*;L*W'i’‘')- 


Using the previous estimates (I5.10p and (I5.12L we deduce the result announced in 
Lemma 15.21 □ 


5.4. Estimates for the memory function K, solution of (15.51) . Unlike stan¬ 
dard memory model^, the estimate of the memory function in the Doi-Edwards model 
is one of the key point in the local existence proof. In fact, the Doi-Edwards model is 
strongly based on the expression of the memory and on the equation that it satisfies. 
We introduce the norm of the memory function that we will controlled in the next 
lemma: 

lll-^lb : = ||-f^||L'»(0,t*;L:»L“ J + ||^Tl^||L<»(0,t*;L”L;j>„) 

-)- ||e^(™''^“*)dTi^||L'»(0,t*;L”L“ J + ll^rVLf||ioo(o_t*;i:^i;,9_^^i|2.«,^) (5.13) 

+ ll<5tItr||i2(o,t*;L*L2 j -I- ||dsiL||LQc.(oy*;L*L2 j. 

^The term “standard” refers to models like K-BKZ in which the memory is a given function, 
usually on exponential type: K{T) = e~^. 
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Lemma 5.3. Let t* > 0 and q ejd, +oo[. 

If there exists /x > 0 such that e^^^'^drKo e and e LyL| g, 

if dgKo e j, satisfies BtLIq ^ 0 , 

ifve L^( 0 , t*-, is free divergence, 

ifdsgeL'^{0,t*;L^J and dsVg e L'^{0,t*-,LU, 

then the problem (EH) admits a unique solution K satisfying drK ^ 0 and 
lll-^^IIU < P3{\\Ssg\\L^{0,t*-,LS^„), l|5sVg||i2(o,t.;L?^^), ll«llL2(o,t.;vy^.<»)), 

where the function depends on q, We, /i and Kq. Moreover this function is con¬ 
tinuous and nondecreasing in each of its variables. 


Proof. The proof is composed of multiple steps since we need to control each term 
of the norm |||iL |||3 introduced by (15.131) . 

•/Step 1: control of iL - Since dsg e L^(0, t*; L“s) and we can directly apply the 
lemma l475l given in the subsection 14.41 We deduce that for almost every {t,T,x,s) e 
(0,t*) X R+ X T X 5 ) we have the estimate 

min{l, inf iLo} / K{t,T,x,s) / max{l, supiLo}- (5-14) 


We remark that the assumption 6 implies a /“-bound on Kq: 


\Ko\\l^^{t)-i / r||dTifol|L-.(r')dr' 
Jo 

’ Jo 


< 






TJ^oIIL^L^, 


(5.15) 


To obtain the first inequality above we used the inequality dr (|| Kq || loo ) / || drKo || l® 
which simply comes from to the inequality sup( 6 ) — sup(a) / sup (6 — a). 

As a consequence of (15.151) . inequality (15.141) effectively gives a /“-bound on K. 


/Step 2: control of BtK - Introducing m = —drK and using the fact that g 
and V do not depend on the variable T, we remark that Equation (15.51) for K is 
translated into a similar equation for m: 

d*™ + + g dgm — djm = 0, 

2Ile 2IIe " 

with the following boundary and initial conditions: 

'ITT'\t = 0 ^ "^L=i ='^L=-i = 

We can apply the lemma 14.51 again to deduce that for almost every {t,T,x,s) e 
(0/*) X R+ X T X (—we have 

0 / inf(— cIt/'o) ^ m{t, T, x, s) ^ sup(—dr/lo)- (5.16) 

Obviously the /“-bound on implies a /“-bound on BtKo so that the 

inequality (15.161) corresponds to a /“-bound on m = —BtK. This estimate (I5.16P 
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also proves that m > 0 . 


•/Step 3: exponential control of BtK - The exponential decreasing is ob¬ 
tained introducing fh = — Cq e'^* where Co = sup {—dTKo(T,x, s) 

T,x,s 

The quantity rh satisfies 

^t'm+-^dTrh-fim +gdsm-= 0, (5-17) 

yJJe iilie 

with the following boundary and initial conditions (we note that we use the assumption 

drKo / 0 ): 

Testing the equation (15.171) with the positive part of to, we easily deduce that to / 0 
(see the proof of Lemma 14.51 for similar result). This implies the desired bound for 
almost every (t,T, a:,s) e (0,t*) x K.+ x T x (—^, 5 ): 

e'" (2H' T, X, s)^ Co. (5.18) 


/ Step 4: control of dr'^K - To have a bound on || Vto||^,^j ;^9 we first note 
that if we introduce to = e™'^^TO, for any A > 0, then we have 


^ + GO ^ + 00 

||VTO|L,^(t,r)dT= e-=»'^^||VTO|U,^(<,T)dr 

Jo ’ Jo 


< 


< 


a +GO . 1 , C + GO . 1 

g-22ITeArjj,^2M ||Vto|| 2 , _^(t,r)dT)' 


||VTO||i,^^(t,T)dr 


It is then sufficient to control ||Vto ||^^^9 The quantity to satisfies 

djTO + ——dyTO — Ato + qdsffi — r-—d^TO = 0, 
aue ^ aue * 


(5.19) 

(5.20) 


with the following boundary and initial conditions: 


it=o 


= —e 


2Hc AT 


StKq, 


'ir=o 


= 0 , 


1 = w _ 1 = 0 . 


Derivating the equation (15.201) with respect to the spatial variable x we obtain 


djVTO + Vn • Vto -f —r—dr'^m — AVm -I- dgfhVq + g dgS/m — -r—d^Vm = 0. (5.21) 
arie ® ^ * 2]r;e s ^ ^ 

We take the inner product of this equation (15.211) by g] Vto|‘?“^Vto, and we integrate 
with respect to a: e T and s e (—^, ^). We deduce 

dt||VTO||«,^^ + ^drllVTOll^,^^ - gA||VTO||^,^^ 

+ ^ f f |VTO|'^“^|dsVTOp + ^ r r |Vto • dsVfhp |Vto|'^“'‘ 

anejTj_i aue JtJ-a 

/ ylo + All ^ 2 , 
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where 


1 . 

Aq = ~q \ (VtJ • Vm) • Vm 

Jt J-i 
1 

Ai = —q ^ J 1 

A 2 =—q j J ^ ^(5sVm • Vm^ 

•/The first term Aq is directly estimate as follows 

Idlol < 9l|VUi|L^||Vm||«,_^. (5.22) 

/The term yli is estimate as follows: we integrate by parts with respect to the 
variable s in order to write Ai = An + An + An with 

^11 = — 9 J J ^ • Vm^ 

^12 = —9 J J ^ 

^13 = —9(9 ~ J J j (^5 ' • Vm^ | VTO|*“‘*m. 

To estimate the term An (which only depends on t and T) we introduce the 
bound on in: 

l^iil < 9l|w||La> f r |dsVg| 

’ JtJ-1 

Due to the Holder inequality we obtain 

Idliil / q\\m\\L^JdsVghijVm\\ll\ (5.23) 

In the same way, the term An is treated as follows: 

|-4i2| / 9 ||w||loo f r |(5sVm| IV5I 

Using the Young inequality and next the Holder inequality, we write 

1^12![ |Vm|«“2|5^Vm|2+^^||m||ig, f f iVg^Vml^-^ 

2WcJrJ.i 2 --JtJ-I 

^ 2^ JtI 1 

Similarly, using the Young inequality and the Holder inequality, the term An is 
controlled as follows 




I Q{Q II II2 ||Y7 ||2 iiYy'^liq'—2 
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Using an integration by parts with respect to the variable s, the contribution A 2 
simply reads 


A 2 = 



‘^dsg. 


Hence we have the estimate 


1 ^ 2 ! < 

/All these estimates (j5.22|l . (I5.23I1 . (I5.24L (15.2511 and (15.261) imply 


(5.26) 


9-1 




llmllij, II Vyll^, IIVm||%^ 


+ '7l|Vv||L*)||Vm||^, 


Multiplying by f ||Vm||^,'^ we also deduce 


dt\\yfh\\li^, + ^^r||Vm||^,^^^2||m||L®J|dsV5||ij_J|VTO||ij_^ 

+ (g-l)2Be||m||i„J|V5||i,^ 
2 


+ -(Il<^s5lli”, + 9A + g||Vi;||L^)||Vm||2,^^. 


We integrate with respect to T e (0, + 00 ). Since mlr^o = 0 the contribution due to 
the term dT||Vm||?<, is non negative. We obtain 


p +00 

^ ’ JO 




||Vm|| 


Ll 


dT 


is? 


r+GO 

+ iq-l)W4Vg\\l._ 

’ JO 

^ + GO 

«, + 9 A +qjlVulli*) J ||Vm||^j dT. 


We now use the bound ||m||L<x,^(t,T) ^ >-)t obtained in the third step 

of the present proof, to deduce 


(?t||Vm|| 


l'LlI 


< 


r + OD 

2Coe^‘||d,V5L., e-™'('^-^)^||Vm|L,^dT 

Jo 

P + 00 

TCo+g-l)2nee2''*||V5||2, ^-2m3t(i.-\)T 

Jo 

p + GO 

+ gA + g||Vu||i^) J ||Vm||2,^^ dT. 


-lITsf/llis?, 













Existence for the Doi-Edwards polymer model 


23 


Using the Cauchy-Schwarz inequality, we obtain for 0 < A < /r the inequality 

^ 2Co e^‘||<5,V5||u,. ( 2 ggef!, - A) ) ' 


UlLl 


2 


Using the Young inequality we deduce 




< 


LlLl 






(5.27) 


We note that ^ ^ Indeed, since 5|s=o = 0 we have 

Vg{t,x,s)= f dsVg{t,x,s')ds'. 

Jo 

By the triangular inequality and next by the Holder inequality we conclude that 

•^“5 

By assumption we know that || Vm||„(0) = IN™'^^<?TVitro||L|,L« „ is bounded for 
0 < A ^ We then now choose A = ^ so that equation (15.271) becomes 




with 


1 2 

a{t) = - + -WdsgU's, + M + 2||Vu||La> 

2 g 

and 6(f) = + 9 - 

We conclude using the Gronwall lemma and the regularity which is assumed for g 
and v: 

^ C'oexp(^ J ^)'^J a(f')exp(^J^ 6^, (5.28) 


with 


C'o = ||e™''^^/2dTVXo||L|LS,.- 

The estimate (j5.28|) shows that Vm G L'^L^ g). From the inequality (I5.19p . 

it implies a bound on Vm in r®(0, f*; s). We also remark that we could have di¬ 

rectly obtained a bound on Vm in T“(0, f*; L\L'^ ^) taking A = 0 starting form (j5.20l) . 
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•/ Step 5: control of dtK and d^K - To get the bounds on dtK and dgK we 
use dtK as test function in Equation (15.51) satisfied by K. We obtain 




J J' (u-Vif)d, 


K 


1 


^Nl ^\\VK\\LiJdtK\\Li^ 

+ \\9\\Ls^J^sK\\Li^J\dtK\\i^2^^ + ^^\\dTK\\L2j\dtK\\L^^^. 

q-2 

Using the fact that L® cz Lj” and using the Young inequality we deduce 


Taking the supremum with respect to the variable T e (0, +oo) we obtairQ 


^ ^MhW^KWl^Li, 


+ Sllffllia) ||dsiir|||^i2 + -z^\StK\\^j^2 

VB.S T aj.s 2 ne ^ 


(5.29) 


By assumptions, we have IIu|||oc. e L^(0, U) and || 5 |||_oo eL^(0,t*) (note that (/|s=o = 0 
so that ||g||L<» ^ II^sSIIl'” )• Using the steps 2 and 4 respectiveljH, we know that 


WhLl^, and llVXlli 
Gronwall lemma to deduce from the inequality (I5.29|) the following bound: 


StKW^^^ and ||VitTH^ooj;^? belong in L®(0,t*). Consequently we can apply the 


ll^i-^llL2(0,t*;L“L2 j + II || ^ 00 ( 0 , 4 *;L*L 2 j < C, 
where C depends on the previous bounds. 


5.5. Proof of Theorem 13.11 For any t* > 0 we introduce the Banach space 

^{t*) =L’'(0,U;1U^’‘?(T)) 

X C([0,U]xi?+;L«(T)) 
xC([0,U]xi?+;L9(Tx(-i,i))) 


^Note that we use the inequality (supj-F(t, T)) ^ sup-j- T)) which simply comes from 

to the inequality sup(6) — sup(a) ^ sup(6 — a). 

^Due to the condition VA'|^_q ^ 0 the step 4 implies that UVA'Il^^g G ^^(0, t*; 
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and for any R* > 0 the subset 

= {(U,G,ir) e 3§{t*) ; 

V e L"(0, dtv e R-Hg), 

dtG, drG e L^iO.R; L^L^), 

K, drK e L”(0, t*; J, dtK e R-L^LIJ, 

dr Vi? e L® (0, t*; ,, n J, d.i? e L”(0, R; L^L^,), 

V|t=o = vq , G|t=o = Go, G|t=o = 

K\t=o = Ko, K\t=o = 1, =0, 

IIHi^ii*, IIIGIII 2 < ii*, |||i?|||3 < i?*}, 

where we recall that the norms ||| • |||i, ||| • |||2 and ||| • |||3 are defined by (lO) . (lOI) 
and (j5.13|) respectively. 

It is important to remark that such a set is non-empty, for instance if R* is large 
enough. More precisely, if 

R* ^ max(Fi(0), |||Go||| 2 , IHi^olb) (5.30) 

then for any t* > 0 we can build a velocity field v* such that (v*. Go, Kq) e JR{R,R*), 
see an example of construction in m p-6]. 

Remark 5.1. If {v,G, K) e JR(R,R*) for some R and R* then the velocity 
field V, the tensor G and the function K are continuous with respect to the time t, the 
age T and the length s. In fact, these continuity properties follow from the Sobolev 
injections of kind bb^’“(0, <*; V) c C([0,t*];X), hold for a > 1. Moreover, they make 
sense of the initial conditions v|t=o = vq, G|t=o = Go and G|t=o = 

More precisely, if (v, G, K) e JR{R,R*) then we have 

v|t=o = Vq in Hq, 

G\t=o = Go m L’fiLl, 

K\t=^ = K^ in 

Noting that L’'(0, t*; L®L^) c: L®(R+; L’’(0, t*; L®)) we also deduce that 

G\t=o = S in M{0,R-,LI), 

K\t=o = 1 m L”( 0 ,f*;L”J. 

Similarly, we have for almost every (t,T) e (0,t*) x R+ the relation K e and 
dgK G L® j,. We deduce that 

KU_^=KU^=0 in L^{Q,R-L^Ll). 


We consider the mapping 

$ : JR{R,R*) 
{v,G,K) 


m{R) 

iv,G,K), 


®This result is based on the following inequality supy ( Jq f(t, T) dt) g: supj. /(f, T) dt for all 

positive function / £ L^(0, t*; Lf). 
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where v is the unique solution of the Stokes problem (15.11) with 

/ = —iHeU • VU + div (T, (5.31) 

where G solves the problem (15.41) depending on v, and where K is given as the solution 
of (15.5p with 


d^g = Vv:S. 


(5.32) 


Lemma 5.4. If {v, G, K) e and r > 2 then we have 


where G is a continuous function vanishing for t* = 0. 


(5.33) 


Proof We will estimate each term of the left hand side of (|5.33|) . 

•/ Step 1: Control of the source term f - The source term f being given 
by the relation (|5.31|) . we directly have 


ll/llL-( 0 ,t*;L«) < • Vi;||Lr-(0,t*;L5) + II div(Ti|ir(o,t*;LJ). 

We follow the ideas of m, generalized by [2] to the d-dimensional case, in order to 
treat the bilinear term v ■ VU. We have 


VU|U.(o.t*;LS) < GP^R 


q + d 
2q 


3q — d 


3q — d 

+ GP^ 


— ” TO* 


i?* 


(5.34) 


The last term div cr is controlled using the orientation tensor S, see (lOD : 


divcrj|i.(o sS a;||VS'||i.(o,t*;L«_^). 


By the definition of S, see (lOl) . the gradient VS reads (making appears m = —dxK) 

^+co _ 

VS{t,x,s)= Vrn{t,T,x, s) ^ y{G{t,T,x)) dT 

Jo 

r+aj _ _ 

-f m{t, T, X, s) y{G{t, T, x)) : VG(t, T, a:) dT. 

Jo 


Jo 

Its g-norm can be estimate as follows 

- + CX) 


IIV5II 


r + aj 

,(i)<^ao ||V?n|U.^(t,T)dr 

Jo 

J ' + OO 

HU»^(t,T)||VG|Uj(t,r) dr. 

0 


(5.35) 


We can read this inequality (15.351) as 
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Note that since ^ R* we have, for all t e 

rOO 

’ Jo 

rOD 

Jo 

^ AiSUe ■ 

We deduce that for all t e (0, t*) we have 

!|V5|U,^(t) + (5.36) 

IWC fi 

Taking the L^’-norm we finally obtain 

II div^|U.(o,t*;L.) < C{t*^R* + (e’-'^** - (5.37) 

where the constant C does not depend on t* nor R*. 

■/ Step 2: Control of the velocity v - We use the following classical result [29]: 
if X c 1" are two Banach spaces with compact injection then the following injection 
is continuous: 


{ve L^{Q,t*-X) ; dtv e ■,¥)] c <^(0, t*; [X, F] i). 


Using X = and Y = we obtain 








'll II 2 

L2(0,t*;L«) 




Vtvp 




(5.38) 




More simply (in fact using the Holder inequality and the continuous Sobolev embed¬ 
ding ^ we have 


ll'’^llLi(0,t*;Wj’«) ^ ^ ll'*^llL’'(0.t*;rV,^''') ^ ^ 


'R\ 


ll'*^llL2(0,t*;Wi'‘”) ^ ll'*^llL'-(0,t*;W,^''') ^ R ■ 


(5.39) 


•/ Step 3: Control of the source term g - Using the definition (15.321) of dsg 
we directly obtain 


\\^sg\\L^{0,t*-,L^„) ^ ll^'*^llL2(o,t*;LM)||‘S'||i<x>(o,t*;L® J. (5.40) 

The velocity contribution is controlled using the Sobolev embedding W^-'^ ^ L®: 


IIVUI 




^ C||n| 




^ Ct* 




^ Ct* 


i?* 


(5.41) 


The contribution of S in the inequality (I5.40p is estimate using its definition (|5.3I) and 
using the inequality I^tXI ^ which comes from the bound |||X |||3 ^ R*: 


ll‘^llL'»(0,t*;L” J ^ C R* e 


fit* 


(5.42) 
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The estimates (I5.41|) and (15.421) allows to write the inequality (15.401) as follows 

(5.43) 

We complete this step considering dgVg. We have 

ll^sVg||i2(0,t.;L9 sS|| VU|| L2(0,tUL*) I V5'||ic»(0, 

+ II ^^'*^11^2(0,4* ;L®) Ill'll L<»(0,t*;L;j'„)- 

We use estimates (15.411) and (|5.42|) . supplemented by estimates (the first comes from 
to estimate (15.361) . the second is a direct consequence of the bound |||l7|||i ^R*)'- 

l|V^tj||L2(0,t.;I,9) ^ ^ Ct~R*, 

to deduce 

^ CR^R*\l + e'^‘*(l + R*)). (5.44) 

Gathering the estimates (15.341) . (15.371) . (15.381) . (I5.39|) . (15.431) and (15.441) we conclude 
the proof of lemma 15.41 □ 

^-Invariant subset - Consequently, using the lemmas 15.1115.2115.31 and 15.41 we 
deduce that if (U, G, K) e R*) then its image {v, G, K) = $(tJ, G, K) satisfies 

lllullli ^Fi(G(G,i?*)), 

|||G|||2<F2(G(G,i?*)), 

|||if|||3<i^3(G(G,i?*),G(G,i?*)). 

To ensure that the ball 3§{t*) is invariant under the action of 4>, we must find t* 
and R* such that 

Fi{G{t*,R*)) sS R*, 

F2{G{t*,R*)) ^ R*, (5.45) 

F3{G{F,R*),G{F,R*)) sS R*. 

Notice that if we choose t* = 0 then the previous inequalities (|5.45|) become (the 
function G vanishes for t* = 0) 

Fi{0) ^ R\ 

F2{0) ^ R*, 

Fa (0,0) ^ R\ 

By continuity argument, taking 

R* = max{2Fi(0), F2(0) + |||Go |||2 , F3(0,0) + HlFolb}, 

we deduce that there exists t* > 0 such that (15.451) holds. For such a choice we have 
the inclusion $(.^(t*)) c 

Note that at the same time , we chose R* so that the inequalities (15.301) hold, this 
insures that J^{t*,R*) 0. Moreover the function $ is continuous and J^{t*,R*) 

is a convex compact subset of ^(F), see [17] for similar properties. We conclude the 
proof using the Schauder’s theorem. □ 
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6. Proof of the uniqueness result. This section is devoted to the proof of the 
theorem l3.2l As usual to prove an uniqueness result we take the difference of the two 
solutions indexed by 1 and 2 such that 

dTVX,eL”(0,t*;LKJ. 


The vector v = vi — V 2 , the scalars p = pi — p 2 , K = Ki — K 2 and the tensor 
G = Gi — G 2 satisfy the following: 


9lc{dtV + Vi ■ Vv + V ■ VV 2 ) + Vp — (1 — uj)Av = div cr, 
div V = 0, 


1. 

r{t,x)=UJ^ S{t,x,s)ds, 


(6.2a) 

(6.2b) 

(6.2c) 

(6.2d) 

(6.2e) 


S{t, s) = ( - StKi {J^{Gi) - S^{G2)) - drK d^{G2))dT, 

dtG + vi ■ VG + V ■ VG 2 + r^drG = Gi ■ Vu + G ■ Vv2, 

2JJe 

dtK + v^-VK + v VK 2 + + (Vui : J 

+ (vui:J + (Vu : J S2)dsK2 = 0. (6.2f) 

together with zero initial and boundary conditions. The uniqueness proof consists in 
demonstrate that u = 0, G = 0 and K = 0. We will initially provide estimates on 
these three quantities. More exactly, we introduce the following quantities 

W{t) = \\Vv\\l,, 

X{t) = 9\c\\v\\l,, 


ni)= Ilmilli.JIGIli.^dT, 

JO 

m = 

JO 


(6.3) 


NIL dr, 


and we will obtained some relations between them. 

Velocity estimate — Taking the inner product of the equation (I6.2al) by v in L 
we obtain 

Die r r 

—dt\\v\\‘j ^2 + (1 - uj)\\Vv\\j ^2 = - cr 'Vv -^c \ {v ' VV 2 ) • v. 

2 ^ ““ JT JT 

From the Cauchy-Schwarz inequality and the Young inequality, we obtain 


Dledt\\v\\l2 + (1 - a;)||Vu|||2 ^ 


1 


1 — CJ 


Ills + 2ine||Vu2||L*||^’||i2. 


(6.4) 
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Introducing ai{t) = 2||Vt)2||L;j> this estimate (16.41) reads using the notations X and W 
introduced below (see (16.dp i: 

X'+ (l-uj)W ^ ^—\\a\\j 2 + aiX. (6.5) 

1 — UJ 

It is important to notice that, due to the assumptions (EH), we have oi e L^(0,t*). 

Stress tensor estimate — From the definition of the stress tensor cr in the 
System (16.21) we have ||crj|i 2 sg w|| 5 ||^ 2 ^ and due to the definition of S (see equa¬ 
tion (|6.2dp h we obtain 

^ + 00 

Ml. \\mMLsJ\nGi)-nG2)\\l. 

Jo 

r+00 

+ M\ \\m\\l. MiG2)\\l^dT, 

Jo 

where rrii = —drKi, i e {1,2} and m = mi — m 2 . By assumption, the function 
and are bounded by and respectively (in practice, recall that we work 
with the function see the preliminary sectionEH- In particular we have M{Gi) — 

y{G 2 )\ < = 5^00 \Gi — G 2 I. In term of Y and Z (see their definition given by (16.31) ') we 
deduce 


klli^ + 

Deformation gradient estimate — Taking the inner product of the equa¬ 
tion (j6.2ep by G in we obtain 

2^t\\G\\l2 + ^ J ^ 

+ r (G-V-U2)-G - r (r;-VG2)-G. 

Jt Jt 

Using the Cauchy-Schwarz and the Holder inequalities, we have the estimate 

2^t\\G\\l2 + 

+ II Vn2||L“ ll^lli^ 

+ ll'nll^^ II VG2||Lrf ||G||l2 . 

2d 

Due to the Sobolev continuous injection ^ and the Young inequality, we 

obtain: 

<5t||G||^2 + ^^^Tl|G'||i2 ^ ll^’llffi + n2(i)||G|||2, 

where the function 

02 ( 1 ) = 2||Gi||^^^|» -I- 2|| Vn2||L^ -I- ll^^2|liMi<i , 


(6.7) 
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and where the constant C depends on q and d. Multiplying this estimate (1121) by 
where mi = —BtKi, and integrating for T e (0,+oo) we obtain 

1 

+ l|mi||i.yTl|G||i.dr^c(x(t) + ir(t))+a2(t)y(t). (6.8) 

■ ■ ■ ■ . ^ 


We note that the constant C contains the value ” ||wi|lia) dT which is bounded 
according to (EH). 

Using an integration by parts we prove that the integral is non-negative. More 
precisely we use the fact that mi\j._^ = 0 and the fact that M -.T ^ is 

a decreasing function. Indeed, we simply use the following property of the supremum: 
for any 0 < T < T' we have 

M{T') — M{T) ^ sup T',x,s) — mi{t, T,x,s)). 

33 ,S 

Since drmo ^ 0 we have seen that dxmi ^ 0 (see the discussion page |S] before the 
statement of the theorem 13.21) . We conclude that M{T') — M{T) ^ 0. 

Consequently the estimate (16.81) now reads 

Y'^CX+ biW+ a2Y. (6.9) 

We also note that bi = C and, due to the assumptions EH, we have 02 e L^{0,t*). 


Memory estimate — We derivate the equation (I6.2fl) with respect to the age T 
and next we take the inner product with m = —drK in Using integrations by 

parts we obtain 



2 

Ll 


+ 


22IIe 


dr II TO 


2 


-f 


QUe' 


sTO 




Lf: 


(Vdi : Si)m^ 


+ 


+ 



m 2 dam 


By the Cauchy-Schwarz inequality we obtain the following estimate: 


II"^IIL 


-f 


22IIe 


drlh 




\\Ll 


< 


llVv 


111^2 




Il"i|li2 


+ l|V^l|k”ll‘S'||L=,J|TO2||L^^J|TO||i2^^ 

+ II II Ll II S 2 IIL* , II TO2 IIL*,, II to|| l2 

+ l|Vt>l||L^||S'||L2 J|TO2||L^,J|d5TO||L2_^ 

+ II Vt’lUj ||5'2 ||l^_^ I|w2||l*, l|dsTO||i2^^ 

+ ll^ll^^ l|VTO2||Ld J|TO||i2^^. 
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From the Young inequality 2ab ^ + b^ and the Sobolev embedding for the last term 

we deduce 


TO 


2 


+ 


SUe 


< 


Czit,T)\\m\\l. 


+ ||Vl)i||Lg>||m2||La> J|S'|||2^^ 

+ II‘^ 2 ||l* J|m 2 ||La> J|Vv ||^2 

+ 2 Uc||V'Uill^^co|| to2||/^oo ||*S'||^2 


+ 21Je||S2||i.J|TO2||i.J|V7;||i2 

+ C'll VTO 2 ||id_ Jib’ll ^ 1 , 


where the function 


( 6 . 10 ) 


Cz{t,T) 




Integrating the estimate (I6.10p for T e (0,T) we deduce (we also use the fact that 


Z'{t) ^ Cz{t)Z{t) + ||Vi;i|U.||m 2 |L^i.J| 5 ||i 2 _^ 

+ 2Ife||Vt!i|||a,||TO2||i2^i^J|5||i2 ( 6 . 11 ) 

+ Wey^\\m2\\li^L^Jm2\\l2^L^W{t) 

+ C\\Vm2\\l.^^,JX{t) + W{t)), 

- -- 2 

where Cz{t) = sup Cz{t,T). Recalling that || 5'||^2 < Y + ^^Z, we can write 

T€E+ _ 

the previous estimate (16.1111 as follows 

Z' ^ asX + ttiY + a^Z + 62 lY, (6.12) 


where we have defined 

a 3 (t) =C'||VTO2||^2^^d^^, 

a4(t) =C{\Xvi\\L-s\\m2\\L\.L^^^ + I|Vr>i||ig,||TO2||i2,i^ J, 

05 ( 1 ) =C(||Vt»i||L£>||TOl||L;^L®. + \\^‘Vi\\L<-£\\m2\\L\.L'S,, 

+ i|Vl’i||L®||m2||L“Lg>„ + I|w2||l;^L;J'J|™ 2||L®Lg',, 

+ l|Vi^illi-ll™ 2 |li^i^^^ + 1 ). 

The assumptions (16.11) imply that 03 , 04 and 05 belong in L^(0,t*). The last contri¬ 
bution makes appear the function 62 which is given by 

h2{t) = y^\\m2\llroo +Wzy^\\m2\ll ra= \\m 2 \l 2 roo + C* II VTO 2 || ^2 . 

_/a3,s j«E,s ja3,s ± »D ^ s 

The assumptions (16.11) imply that 62 e L“(0,f*). 
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Uniqueness result — Finally, we perform the following combination: 

(1 - w)(l631) + (EH) + eEH + eEmi) 

using e > 0 small enough to control the quantity W{t) of the right side member by 
the same quantity on the left side member. This is possible since the functions bi 
and 62 are bounded. More precisely, we use 

1 + 61 + ||^2||L®(0.t*) 


We deduce an estimate of kind 

(X + Y+Zy £ia{C + Y + Z), 

where the function a is a linear combination of a^, i e {1,2,3,4}, and consequently it 
belongs in The classical Gronwall lemma and the initial conditions X(0) = 

U(0) = Z{0) = 0 imply that X = Y = Z = 0. We deduce that d = 0, G = 0 
and K = 0, that concludes the proof. □ 


7. Global existence. In this section, we prove the main result of the present 
paper, that is the global existence of a solution to the System (12.81) in the two di¬ 
mensional case, see Theorem 13.31 The proof is based on the following fundamental 
remark: the stress given by the relation (I2.8cl) is “naturally” bounded. This bound 
implies a bound on the velocity field, from which we can deduce a new bound on the 
stress gradient. The conclusion of the proof consists in to insure that this process is 
consistent. 

The main difficulties are the evaluation of the nonlinear terms that can blow the esti¬ 
mates obtained in the proof of local existence. These nonlinearities are present many 
times in the Doi-Edwards model: 

1. The nonlinearities of the Navier-Stokes equations (I2.8ap - (l2.8bl) are always 
present in the complete Doi-Edwards model. On one hand there is the con¬ 
vection term v ■ Xv which is discussed in many articles. We handle this 
first difficulty using the so-called L^L'^ estimates of the heat kernel. On the 
other hand there is the nonlinear coupling due to the presence of the stress 
term div cr. 

2. Other nonlinearities are present in the evolution equation (I2.8el) for the defor¬ 
mation tensor G. The most restrictive one corresponds to the term G ■ Xv. 
Since we will have relatively little information on the velocity gradient, we will 
prefer to work with the quotient making appear the product \G\y{G). 
The key point is the fact that this product is well behaved even for large value 
of G. 

3. Finally, there are also several nonlinearities in the equation (I2.8fl) describing 

the evolution of memory K. The most restrictive one comes from the product 
{Xv : Jp S)dsK. To solve this problem, we introduce a suitable function ^ in 
order to obtain a better estimate on and therefore a better estimate 

on BtK. 
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Let r>2, (7>2, 7>0 and ^ > 0. We consider the initial data (?;o, Gq, Kq) satisfying 
^^0 6 

drCoeL^Ll, detGo^y, 
e^^dri^o e e^^/^drVKo e „ d,Ko e 

and drKo ^ 0. 


It is then possible to use the Theorem 13.11 and then consider the solution {u,G,K) 
to System (12.81) in [0,t*], which satisfies the initial conditions (12.61) . This solution 
possesses at least the following regularity: 


K, drK, e>^ drK e L^{0,t*;L^L'^,), 

dTViLeL”(0,t*;L^L«^,nL|L«_J, 


dtveL-{0,t*-,Lg), 

dsG, dtGe L’'(0,t*;L“L«), 

dtK e L\0,t*-L^Ll ), 


In the following and as previously we recall that we denote by C constants that may 
depend on the initial conditions, on the physical parameters, on the integers r, q, on 
the bounds oS^oo and 5^'^, and on the time t*. Note that these constants will always 
be bounded for bounded t*. 

7.1. Maximum principle for the kernel K and their derivatives. First 
we resume the proof of the local existence theorem to deduce the following bounds 
(see the estimates (I5.14L (15.161) and (|5.18|) b 

Lemma 7.1. There exists a constant C such that 


\K\ < G, 

0 ^ -drK ^ G, 

0 ^ ^ G. 


7.2. Additional bounds for the stress tensor and for the velocity field. 

The first Proposition that we introduce is one of the fundamental points which give 
the global existence result. In particular, this kind of result is not proved for other 
viscoelastic systems like the Oldroyd models. 

Proposition 7.2. We have the following U^-bound: 

l|l^llLa>(0,t*;L“) ^ G. 


Proof. Recall that the stress is given by the integral relation (see ( 12 . 80 ^ 1 : 

1 

(T{t,x) = uj S{t,x,s)(is. 

J-h 

To have a bound on cr it suffice to have a bound on S which is given by (see (I2.8d() h 

p+oo 

S{t, x,s) = - drKit, T, x, s) y{G{t, T, x)) dT. 

Jo 
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Since the function is bounded on £(R‘^)\{Q} (see Proposition 14.11) . and since the 
function G has values in r(IR‘^)\{Q} (see Lemma we deduce that the composed 
function ^{G) is bounded. 

Moreover from Lemma l7.ll we know that —dTK > 0 and we deduce that 

p+oo 

|5(t,a:,s)| < dTK{t,T,x,s)dT = lim K). (7.1) 

Jo T—>+00 

Using the fact that K is bounded (see Lemma 17.11 again) we conclude that S is 
bounded. □ 

By virtue of this proposition 1 7. 21 we can deduce that the velocity field v satisfying 
the Navier-Stokes equations (I2.8al) - (l2.8bl) have more regular. More precisely we prove 
the following Proposition (the proof is detailed in [3] and i)- 

Proposition 7.3. If the integers r and q satisfy | + y < 5 then for all t e (0, t*) 
we have 


I^ C'ln(e + || Vcr||)), 

II < C'(l + l|V<T||j;^r(o^t;i«)). 


7.3. Control of the spatial gradient of the orientation tensor S. The 

goal is now to use the previous bounds (on the velocity field) to obtain a bound on 
the stress gradient Vcr. In fact, the stress gradient is directly linked to the orientation 
gradient VS (see Equation (j2.8c|) l and we are then interested in this quantity: 

Proposition 7.4. For any convex function f ^ R+ we have: for all 

te( 0 ,U), 

/ - 2 WtuT r N 

where ya and yb are defined by 

m +00 r 

\\VM\\l,Jt',T)dTydt\ 


t r+cc 


Vbit) = f f 

Jo Jo 




VG 


|G| 


{t',T) dTdt'. 


The quantity M is defined by M = f{Th) and m = 


(7.2) 


Proof As previously we denote m = —BtK■ We derivate the stress tensor given 
by the relations (|2.8c|) and (|2.8dl) with respect to the spatial variable: 

y+c» 

VS(t,x,s)= Vm(t,T,x, s) ® y'{G{t,T,x)) dT 

Jo 

'-V-" 

(VS)“(i,a 5 ,s) 

C + GO 

+ 7n{t, T, 0 ^, s) T, x)) : VG(t, T, x) dT. 

Jo 

- - - - ^ 
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We have the estimate 

(7.3) 

so that we will independently control ||(V5')“||^r(o_(.2,5 and ||(VS)^||j;,r(Q 

■/ We use the definition of M = ^{fh), fh = to write the quantity (VS')“ 

as follows: 

r+co p— QUe/iT 

{VSnt, x,s)= --VM(t, T, a:, s) ® y{G{t, T, x)) dT. 

Jo ^'(m{t,T,x, sj) 

We deduce that 1|(VS')“||2,3 ^ is bounded using the L^-bound on y (see the Propo¬ 
sition ED), the bound on (by the convexity assumption and the positiveness of fh 
we have : ^'{fh) > C^(0)) and the triangular inequality: 


:|(V5)“L, (t)<^o 


I 


+ 00 ^-WeiiT 


Jo e'(o) 

We next use the Cauchy-Schwarz inequality and obtain 


VM||i,,(t,r) dT. 


a 'OO -2WtllT 1 , r + GO ^ 1 

. TTT") (Jo ■ 


Jo m 

Taking the L’'-norm on (0, t) we deduce 


ll(V5)“irLqo+;L£,.)<=^»l 


a 'OO -223Jc/ir s 

-dT) 

0 f'(0)2 / 


dT y,(t). 


Jo e(oy 

/ Due to the Proposition 14.II we control the last contribution (VS)^ by 

VG 


(7.4) 


(VS)'’(t, X, s) ^ =5^^ m{t, T, X, s) 

JO 


0 

By the triangular inequality we deduce that 

r' + OO 

Jo 


IGI 


(t, T, x) dT. 


r + oo 

l(vs)^IU^_(t)^y^ l|m|U^^(t,T) 

Jo 


VG 


|G| 


(t, T) dT. 


From Lemma [7.11 there exists a constant C such that ^ Ge^d WcT)^ 

We deduce 


^ + 00 

ll(VS)'Lj,(t)<G . 

Jo 


Mt-WcT) 


VG 


|G| 


{t, T) dT. 


Taking the T’'-norm and using the triangular inequality, we obtain 






VG 


|G| 


\ !/’■ 

(f',T)dt') dT. 
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Writing e = e iiWcT], using the Holder inequality, we obtain 

\ r—1 




C + 00 

Jo ' 


-/x2BeT 


n + CO rt 

JQ Jn 


gi'/it'e-M^DcT 


VG 


|G| 


(t',r) dt'dT. 


This inequality (17.31) corresponds to 


(7.5) 


/ Finally, we add the contributions (17.41) and (17.51) and we use the inequality (17.3p 
to conclude the proof of the Proposition 17.41 □ 


7.4. Control of ya- The goal is now to analyze the quantities ya and y^- This 
subsection is devoted to the control of y^, the next subsection will be devoted to the 
control of yi,'. 

Lemma 7.5. There is a family of functions f, parameterized by k ^ 1, for which 
the function ya introduced in the Proposition \7.4\ satisfies 

y'a < <^(1 + ln(e +y)ya) + + \n{e + yf yr)yl~^. (7.6) 

The function y is defined by y = k~ 2 ya + yj, and the constant C does not depend 
on fc > 1. 


Proof. This proof is decomposed into two main steps. In the first step we establish 
the evolution equation satisfied by the quantity VM which appear in the definition 
of ya- In the second step we show that for a ’’good” choice of the function f we 
can deduce an interesting estimate for VM, that corresponds to the estimate (TLel) 
announced in the lemma FfTSl 

/Step 1: evolution equation for VM - We recall that m satisfies 


+ gdsm - ^d'^m = 0 , 

iLUc iLUc 

with the following boundary and initial conditions: 


(7.7) 


and where the function g is given by 

g{t,x,s) = Vv{t,x) : f S{t,x, s') ds'. 

Jo 

To obtain the equation satisfied by ^(m) we multiply (17.71) by We deduce 

dt(C(m)) + ^dT(C(m)) + gds{i{rn)) - = 0. (7.8) 

One of the key points of the proof is that we want the function f provides more 
regularity. After multiplication by f,'{rn) the regularizing term d^m has now becomes 


r'(m) 

(f (to))2 




^'(m)d^TO = (9f(^(m)) 








38 


LAURENT CHUPIN 


Hence the equation (EH) writes 

+gds{am)) - ^ = 0 . 

(7.9) 

Multiplying this equation EH by we obtain the following equation on the 

quantity M = e™'^^^(m): 


d.M + 

Finally, denoting m = e™'^^TO and ■C(w) = e®^'^^^(m), we obtain 


(d,M)" = 0 . 


dtM + -^M + gd,M- = 0. (7.10) 


QUe 


Taking the spatial gradient of (17.101) we obtain 

djVM + Vv.VM + — uVM + VodgM + gdsVM — 

QUe y s y s 2jje « 


1 


^"(m) 


VM ~^2 . 2 ^"(m) 


(7.11) 


■(d.M) + 


■dsMdsS/M = 0 . 


'?'(™) 2 ne (^'(77i))2 

Remark that the first line of the equation (17.111) is exactly the same that those obtained 
during the local proof (see equation (15.211) '). The introduction of the function ^ makes 
appear the two last terms in E-HP and at this stage we hope that such terms brings 
more estimates. 

/Step 2: choice for the function ^ in order to estimate VM ~ We hrst 
impose that the function ^ is the solution of the following Cauchy problem: 

U”(^) = eo(2/(x) + 6)"(yW 
U(o) = o, 2/'(o) = a. 


(7.12) 


The parameters ^i, C 2 and k are assumed to be positive and will be judiciously 
chosen later. We simply note that these constants do not depend on x and s but may 
depend on the times t and T. Some results about this ordinary differential equation 
are given in the last preliminary, see Section [4.51 
With this choice, the equation (17.111) becomes 

dt VM + Vu.VM + v^drVM — /rVM + S/g dgM + g ds VM — VM 

2ifc 2iJc ^ 

+ y^(C(m) +6)"“'(<5 «M)Vm+ ^(C(m) +6)"5«Md«VM = 0. 

We now proceeding as in the proof of the local existence theorem (see page [ 201 ). 
We take inner product of this equation (17.131) by g|VM|'^“^VM, and integrate with 
respect to a; e T and s e (— 5 , 5 ). We deduce 

1 


dt||VM|| 


LI,. + 


« -qg,\\VM\\l, 


^ ir. l^"r=l<>.VMP + JJ(_ |VM . .l.VMP IVMr* 

+ 1 (^(”^) + '^2)^ ^ I VM|'^|dsMp / Ho + Hi + H2 + H3, 
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where 


^0 = 
Ai = 
A 2 = 
A3 = 


1 . 

q\ I " (Vv • VM) • 

JtJ-I 

1 

9 J J' (vg- 

1 

1 

J J \ (C(w) + 6)'" ( 5 .VM • vm) 


JT J- 
2g$o 

21Je JTJ- 


The quantities Aq and A 2 are estimate by (see (15.221) and (15.261) respectively for the 
same kind of estimates): 


l^lol <9||Vu||ij||VM||«,^, 


l^2| < Wdsghs,^ l|VM||«, 


The proof that we use to obtain global estimate fundamentally differs from the proof 
presented for the local existence especially in the control of the term Ai. We control 
the contribution Ai without integration by parts: we simply use the Young inequality 
and then the Cauchy-Schwarz inequality: 


1 

l^il < + 

■ ■ 1 _ 

JtJ -1 l^ffPlVMl 


qWe 


^(m) + C 2 


q-2 


(7.15) 


Using the fact that the solution ^ of the Cauchy problem (17.121) is increasing (see 
Section 14.51 for an explicit expression of the function ^), and the fact that m is non¬ 
negative (and vanishes), we directly estimate 


1 


1 

^(m) -t 6 


l'(0) + 6 


Using the Holder inequality, the estimate (17.151) implies 


1^1 




22ne 




m) + ^ 2 f-^\VM\'J\dsM\^ 


+ 


T J —i 

m II 

-^11^511 


2 k 


||VM| 


q-2 
Ll ■ 


(7.16) 


The term A 3 is new with respect to the estimate introduced in the local proof theorem. 
This term comes from to the function ^ introduced here. We must then control this 
term with the supplementary contribution given by the function ^ itself. In practice. 
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the term is treated similarly as the term Ai^ that is using the Young inequality 
and then the Cauchy-Schwarz inequality: 


l^3| < 


+ 


2 gCo 

kT5z' 


Iie(m) + C2|lir |VMr-2|d.VM|2. 


Consequently to control A^ with the left hand side member of (17.141) we want to 
choose the parameters ■fi) C 2 and k such that 


2 gCo 

kWe 


Iie(m) + 6lli£ 



(7.17) 


By lemma FTTl we know that for all {t,T,x,s) e (0,<*) x R+ x T x (— 5 , 5 ) we have 
0 ^ m{t,T,x,s) ^ rhoo = C. Since the function ^ is increasing the condition (17.171) 
also reads 


Using the preliminary result given by the Proposition l4.61 we have an explicit expres¬ 
sion for the solution making appear a function F. The condition (17.171) is then 
equivalent to 

In the sequel, we choose ^0 and with respect to the parameters ^2 and k as follows: 


With this choice, the inequality (17.181) holds, hence the inequality (I7.17p holds too. 


Remark 7.1. The choices of and are fundamental since they ensure the 
validity of (17.171) . but they are also for the following three reasons: 

1. First it is important to notice that with such coefficient the real fijh) is defined 
for any rh e [0, TOqo] . Indeed the set of definition of the function f given in 
the proposition \4-t>\ writes 



c pfc + l 

50 S 2 

F+i 



^ ^ 1 

F\ 

[k\ 

if 

Fb-b) 

"‘'00 1 

1 


[ 0 ,moo], 


k ^ 

since we have ^(^( 2 ) ^ ^2 — ^ 2 ^ < i = F{X). 

2. We also note that with this ehoice the estimate (17.161) for Ai becomes 


1^1 




22IIe 






2 k 


V-lIlL 




(7.20) 


We will see later that the coefficient k at the denominator of the last term 
will make this contribution as small as desired (letting k to +co). 
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3. Using the expression of the function F given by (14.81) . the value of^i introduce 
by (IM is written 


1 1 r 

=-e(fc+i)fcfc+i 

Woo Jo 




da;. 


Performing the following change of variable z = in the integral, and 

using the expression (I7.19|) for we obtain 


1 1 r 

Woo Jl 


e (fe+i)fcfc+i dz. 


We can find a lower bound for the integral using the faet that for z ^ fc(|) '=+1 
we have 

_ 1 _ k 1 

0 (fc + ljfcfc + l ^ 0 2(fc+l) ^0 2 


We obtain 


fi > .^e^2- l)e“5. 

Tri(yQ V 2 ) 

Finally, since > 1 and (|)’^ = 5 ; 1 for k ^ 2, we 

deduce 




(7.21) 


where the constant C does not depend on ^2 nor on k. 

Using (17.4L (17.4[) and (17.201) . the estimate (17.141) now write 

+ ^dT||VM||«,^^ <( 9 ||V^|U^ + ||d,5||L-, + qu)WM\\ 

g2rreC|ii^__ii2 ||^~||5_2 


+ 


2 k 


iivffiii, iivMir 


(7.22) 


We multiply this estimate (17.221) by and deduce 


Wzfl 




We integrate with respect to T e (0, + 00 ). Since m| 7 ’=o = 0 the contribution due to 
the term dT||VMB<, is non negative. Finally, we obtain 


( C + 00 \ 2 c + co 

J ||VM||2, ^ dTj s :(2 IIVi;IIL® + -IdspU-s^^ + 2 p) J ||VM||^, ^ dT 

r+® 


+ 


J 

Jo 


(7.23) 


<iT]\\Vg\\U- 
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We finally chose k independent of the time T, and 

1 


6 = 


1 + T’ 


(7.24) 


so that the estimate (17.231) becomes 

dT^ <(2 ||V«||l- + 


q 


L'S. + 


+ 00 


IIVMiii, dr 


2ne|_ .,2 


To make appear an equation of evolution on ya we multiply the estimate (17.231) by 


(C”l|VM|L,^dT)5 \ We deduce 


QUe, 


y'a + -||<5sffl|L*, +ry)ya+ 

Integrating the Equation (17.251) with respect to time t' e {0,t), we obtain 

T 

y'ait) ~ y'aid) <(?’||Vt;||ia)(o,t;L«)) + - || dsff || L>»(0,t;L” J + fy)ya{t) 


(7.25) 




The last integral is treated using the Holder inequality again (we also use the fact 
that ^ 0). We deduce 


V 

y'ait) - 2/1(0) sS(?’l|Vi>||Loo(0.t;L*) + -||dsg||Lao(o,t;L«j + ry)ya{t) 

mr (7-26) 

tWZ ||9 / \ 1 _ 2 

+ ^l|V5llLr(o.t;LJ,,)2/a(<) 

Using the definition oi g = Vv : S ds we note that 

ll'^2/llL>-(0,t;Lj,s) ^ ^II'^’*^IIl‘=°( 0+;L“) ll^‘^llL>-(0,t;Lj,a) + 2||V I’ll ir(Q ) || •S'|| 

(7.27) 

Moreover, using successively the Droposition l7.41 the remark [7.11 (more precisely equa¬ 
tion (17.21^ 1 and the expression (I7.24p for ^2, we have 

W oo -2Wefj.T . \ 

^ -^^dTjyait)^+y,itfrJ 

< C{^ya{t)- + yb{t)-). 

Introducing y = k~iya + yb and using the inequality + bi ^ {a 4- 6)^, holds 
for r > 2, a > 0 and 6 > 0, we obtain 
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We also deduce from the Proposition 17.31 land using (17.281) 1: 

ll^^'*^llL'-(0,t;LJ) ^ C'(l + IIVcr||£,r(o_j.j;,9)) 

^C(l + ||V5|U.(o.t;i.,,)) (7.29) 

< (7(1 + 

From the Proposition 17.31 again we hav£0 

II V'*^llL‘»(o,t;Lg’) ^ (71n(e+ II Vcr||j;,,.(Q 

^ Cln(e + (7.30) 

^ C'ln(e + j/(t)). 

We have already seen that the orientation tensor is bounded (see Equation (ED 
obtained to get a L°°-bound on the stress cr): 


ll'S^llL‘”(0.t;Lg’J ^ (7- 


(7.31) 


Finally, the initial value 2/(i(0) is estimate using the assumption e^^/^Vwo e s'- 


, 1 ( 0 ) = 
*(I 


II Vtoi 


oil Li 


dT 


)' 


r ^ + CO : 

-^^dT)'(J e'^'^llVmoll^, ^dr) 


(7.32) 




Gathering the six previous estimates (|7.27I1 - (I7.32|1 . we write the inequality (17.261) as 
follows 


C , 


2 1-F 


y'a < G(l + In(e + J/) Va) + — (l + ln(e + j/)2 ?/-)y, 


7.5. Control of yt,- We will now control the quantity yb in order to obtain the 
following result: 

Lemma 7.6. The function yb introduced in the Proposition \l.Ji\ satisfies 

y'b < (7(l + ln(e + j/) yb), 

where we recall that the function y is defined by y = k~^ya + yb and the constant C 
does not depend on k ^ 1. 

Proof Recall that the function yb is given by (see (17.21) 1 


n +oo 

) 


VG 


{t',T) dTdt'. 
Li 


^We also use the fact that there exists two constants C and C (depending on r) such that 
for any y > 0 we have ln(e + Cy~) ^ C'Mn(e + y). This is due to the fact that the function 




ln(e + CyT) 
ln(e + y) 


is bounded on (in fact it is continuous and satisfies £(0) = 1, lim£ = —). 

+00 r 
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To prove the lemma 17^ we must write an equation of evolution for The equation 
satisfied by G reads 


VG = G Vv, (7.33) 

where we recall that V corresponds to the operator V = dt + dx- We take the inner 
product of Equation (17.331) by —q\VG\‘^\G\~‘^~'^G-. 

\VG\‘iV\G\-‘^ = -g|VG|«|Gr«-2(G • Vi;) : G. 

Using the Cauchy-Schwarz inequality, we deduce 

\VG\'^V\G\-'^ ^ g|VG|«|Gr‘?|Vi;|. (7.34) 

Next we take the spatial derivative of Equation (17.331) . We obtain the following 3- 
tensor equation 


VVG = VG • Vi; + (G • V^i;)^ - Vi; • VG. (7.35) 

More precisely, the component {i,j, k) of this equation reads 


dydiG jk — diG j^d + G j^d^diV diV^d^Gj^. 

Taking the inner product of this equation (I7.35P by (7|G|“'^| VG|‘?“^VG and using the 
Cauchy-Schwarz inequality we deduce 


|Gr«I?|VG|« ^ 2(7|VG|«|Gr‘?|Vi;| + q\VG\‘i-^\G\-^'^-^'^\V^v\. 
Adding this inequality (j7.36p with inequality (17.341) we deduce 

I?(|VG|'?|Gr«) 3g|VG|«|Gr‘?|Vi;| +g|VG|«-i|Gr('?-^)|V2i;|. 

Integrating with respect to the spatial variable we obtain 


VG 


|G| 


-l- di 


VG 


|G| 


^ 3(7 


I 

Jt 


VG 


|G| 


|Vi;| -I- q 


J 

JT 


IVGI 


We now use the Holder inequality to write 


VG 


|G| 


+ dx 


VG 


|G| 


We multiply (17.371) by 


U 

VG 


< 3g 


VG 


|G| 


I|Vi;||l2> + q 


IVGI 


|G| 


|G| 


9-1 

I 

J-<!E 


9-1, 


|V^i;|. 


V^i;! 


VG 


|G| 


+ dx 


G 

VG 


r—q 


to have 


|G| 


^ 3r 


VG 


|G| 


IIVvIIloj -I- r 


VG 


|G| 


r —1 


llV^i;! 


Using the Young inequality we obtain 


VG 


|G| 


-f dx 


VG 


|G| 


^ 3r 


VG 


|G| 


II Vi; I 




VG 


|G| 


+ (r - 1) 


r-l|lv2i;| 


LI' 


iJM) 


(7.37) 
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,11 VGr 

The quantity Q{t,T) = 9 satisfies 

d^g + drG ^ 3r||Vt;||i^^? + (1 + nr)g + (r - 

We multiply by and integrate for T e (0, + 00 ). Using the property = S 

we deduce that the second term is non negative (it equals ny'Y- We obtain 

y'b ^ 3r llVvIli^y^ + (1 + yr)y[ + (r - l)’-ie'"’'‘|| . 

We can write for any 0 < t' < t the relation 

y'bit') ^ 3r t|Vt>||ia.(o,t;L»)y6(t') + (1 + yr)y'Yt') + (r - l)’-ie'"’'*||(t')- 
Integrating now with respect to time t' e (0, t) we deduce 

y'bi.'^) ~ 2/b(0) ^ C'dVi;||iQo(o,t;L;j>)yh(^) + 1 + yb{t) + (7-38) 

The value of y(,(0) is given with respect to the initial condition Gq: 

VGo 


f 

ym = 

Jo 


^-Wc/J.T 


IGn 


AT)dT. 


We will note that 2/(,(0) bounded since Go 6 L®(K+; and |Go| > 7: 

^ 7^ 2Ilc y ^ 

Moreover, using the relation (17.291) and (17.301) again, that is 


WML^OALD^Cil + y^^Y, 
l|Vv||L=»(o.t;L*) < Gln(e + y), 

the inequality (j7.38|) becomes (note that (1 + ^ 2’'“^(1 + y)) 

y'b ^ G(l + ln(e + 2/) yb), 

that concludes the proof of the lemma 17.61 □ 

7.6. End of the proof for the global result. Adding the results obtained in 
Lemmas 173] and 1731 we deduce that y = k ^ya + yb satisfies 

y' ^ G(l + ln(e + y)y) + ^(l + + yf y-)y^~-. (7.39) 

Using the Young inequality we have y^~i ^ G(1 + y) ^ G(1 + ln(e + yYy). The 
equation (I7.39|) implies 

y' ^ G(l + ln(e + 2/) 2/) + ^ (l + ln(e + yfy). (7.40) 

Since y > 0 we have 1 ^ (e + y) ln(e + y) ^ (e + y) ln(e + yY and the equation (I7.40p 
implies 
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Dividing by (e + y) ln(e + we get 

w'^ C (—TT — w) where w = -r. (7.41) 

Vfc 2 J ln(e + y) 

This linear equation (17.411) . combining with the initial value w{0) = is equiv¬ 

alent to 




4 - 


^-Ct 


Tn(e + yo) 


(7.42) 


We note that the right hand side member in (17.421) can vanish at time to given by 


to = ^ In (1 -f 




e + yo 


It verifies lim to = -l-oo. 

A;^ + oo 

But by construction, we have re < 0 so that the time of existence of the solution y 
satisfies > to- Since the constant k can be chosen arbitrarily large, this time can 
be choose as large as possible: we deduce that y is bounded up to time t* if k is large 
enough. 

This bound on y implies the bounds on S {via (17.281) and (17.311) 1 and the bounds 
on Vt), using (|7.29l) and (j7.30|) . That completes the proof. □ 

8. Conclusion and open problems. In this paper we proved that the famous 
Doi-Edwards model describing the dynamics of flexible polymers in melts and concen¬ 
trated solutions is mathematically well-posed. In particular, we show that the solution 
exists for all time in 211, irrespective of the data (not necessarily small). Such results 
naturally bring a lot of new open questions: 

/ The first one is about the periodic assumption: the results that are proved in the 
present paper correspond to the case where the spatial domain is periodic. The reason 
of this choice is purely mathematics and it is certainly possible to extend some of the 
results to the bounded domain case using Dirichlet condition for the velocity field (we 
note that such a result - local in time - have been obtained for polymer model with 
integral law, see 0). 

/ The second point which can be a source of interest is the existence of stationary 
solution. From another point of view, the theorem 13.31 indicates that the solution 
exists on (0,t*) for any time t* > 0. But we don’t know if there exists a solution 
bounded on ( 0 , -l-oo), whose the long time limit would be a stationary solution. 

•/ Finally, as it was specified in the introduction, much progress has been made on the 
modeling of both linear and branched polymers since the pioneering works of M. Doi 
and S.F. Edwards. There exists many models, usually derived from the Doi-Edwards 
model, which take into account more and more complex phenomena. An example of 
such a model is the so-called pom-pom model introduced by McLeish and Larson [55] , 
and later modified by Blackwell et al. [1]. To the best of my knowledge, no mathe¬ 
matical results exists for such problem and a great challenge would be to show that 
they are globally well-posed. 
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